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ABSTRACT 

We present OASIS observations obtained at the Canada-France-Hawaii Telescope for the 
SB(rs)c galaxy NGC 4900. About 800 specti-a in the wavelength range 4700-5500 A and 6270- 
7000 A have been collected with a spatial resolution of ^ 50 pc. This galaxy is part of a sample 
to study the stellar populations and their history in the central region of galaxies. In this paper, 
we present our iterative technique developed to describe consistently the different stellar com- 
ponents seen through emission and absorption lines. In NGC 4900 we find many young bursts 
of star formation distributed along the galaxy large scale bar on each side of the nucleus. They 
represent nearly 40 per cent of the actual stellar mass in the field of view. The age for these 
bursts ranges from 5.5 to 8Myr with a metallicity near and above 2Zq. The extinction map 
gives E(B—V) values from 0. 19±0.01 near the youngest bursts to 0.62±0.06 in a dusty inter- 
nal bar perpendicular to the large scale bar. The Mg2 and Fe I absorption lines indicate the 
superposition of a background stellar population with an age between lOOMyr to 3 Gyr and a 
sub-solar metallicity on average. We propose that all these episodes of star formation are the 
consequence of a secular evolution. In this scenario, the galactic large scale bar plays an im- 
portant role with respect to the recent bursts and the dusty nuclear bar observed. The iterative 
technique allows us to improve the determination of the stellar population parameters, mainly 
an older age is obtained for the old component and more reliable stellar population masses 
are found. A composite/transition type activity in the galaxy nucleus is also revealed with this 
technique. 
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1 INTRODUCTION 

The central region of spiral galaxies hides a large variety of ac- 
tivities and stellar populations. Observations of nearby galaxies 
and theoretical models for the formation of galaxies indicate the 
superposition of disk stars, old and metal-poor objects from the 
bulge, and often additional populations related to peculiar phe- 
nomenon near the galaxy center (see for example the study o f 
the central region of the Milky Way by Ivan Loon et alj l2003h . 
The assemblage and evolution of these components is still puz- 
zling. The presence of a large scale bar in the galaxy is believed 
to be an efficie nt way to drive gas toward the central kiloparsecs 
of the galaxy jNorman et al.|[l993) , which may trigg er star for- 
mation or turn on the central engine jHaan et al.l l2009). However, 
iRegan (STT euben (2004) showed that a bar alone will only bring 
gas to the central region, but further structures and perturbations 
such as nuclear bars, rings, or spirals are needed to drive the gas 
deeper into the galaxy. Nucle ar rings are found in about 20 per 
cent of spirals l lKnaDenll2003) and these are very often, if not al- 



ways, barred spirals ( lButa&Combeslll996h . Interactions, mergers 
or cosmic filaments are considered t o play a role in the accretion 
of gas and its flow process as well (Combes 2008). Secular evo- 
lution is proposed to be important for the formation of the bulge, 
called a pseudobulge, for late-type spirals ( IKormendv & KennicuttI 
l2004h . This slow process may be internal - due to non-axisymetry 
from a bar, oval disk, global spirals, or a three axial dark halo ~ 
or environmental, due to galaxy harassment for example. It rear- 
ranges the material from the disk, creating inner and outer rings 
and a dense concentration of gas in the galaxy center which may 
in turn form new stars, nuclear structures, and becomes the pseu- 
dobulge. As opposed to true (quickly-formed) bulges, which are 
believed to ass emble in less than 1 Gyr early in the galaxy history 
(^10 Gvr ago; iRenzinill 19991 : IZoccali et alj[2006l : ICunha & SmitH 
1 20061 : iBallero et alj|2008[) , pseudobulges build up on time scale of 
1 to 3 Gyr, are rapidly rotating, and have luminosity profile closer 
to a disk with a Sersic inde x between 1 and 2 (for a summary of 
pseudobulges properties, see I Kormendvll20()3) . 

One way to deepen our understanding of galaxy formation 
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and evolution is to obtain a detailed description of their stellar and 
gaseous contents. In distant galaxies, where individual stars can- 
not be resolved, information about the stellar populations (i.e. their 
age, metallicity, initial mass function [IMF], stellar formation rate, 
mass, etc.) is gathered using evolutionary synthesis techniques on 
the integrated galaxy spectrum. This type of study is somewhat 
difficult and limited. For instance, it is highly dependent on the 
wavelength range studied since specific generations of stars will 
dominate the integrated flux while the emission from others will 
be washed out. For example, OB stars, even in small number, may 
acco unt for most of the UV fl ux and hide older generations of stars 
(e.g. lPellerin & Robertll2007h . In the optical, young and old popu- 
lations may contribute to the spectrum, but often one is ignored for 
a quick analysis of the dominant component. Furthermore when 
considering the addition of several stellar populations, a unique so- 
lution is difficult to obtain depending on the number of stellar in- 
dicators available. This is mainly caused by degeneracy effects be- 
tween the age, metallicity, a nd reddening, as ca n be seen in models 
for a single population (e.g. ICarter et al.ll2009l) . Nevertheless, new 
instrumentation for the observation of distant galaxies encourage 
us to pursue the development of detailed synthesis techniques for 
the reconstruction of the history of galaxies. In particular, integral 
field units, with high spatial resolution spectroscopic capabilities, 
are very promising if one wants to isolate and characterize individ- 
ual stellar populations in crowded regions of distant galaxies. 

In this paper, we present a study of the stellar populations 
and gas in the central region of the barred galaxy NGC 4900 (RA: 
I3h00m39.Is, DEC: +02°30m05s, J2000) using visible data ob- 
tained with the integral field unit OASIS (Optically Adaptive Sys- 
tem for Imaging Spectroscopy). NGC 4900 is a SB(rs)c galaxy lo- 
cated in the Virgo cluster at a distance of 13.3 Mpc (using a red- 
shift of 960±3kms"^ from NED and Ho=72kms~^ Mpc"^). It 
is located at about 12° South-East of M87. It does not show any 
evidence of a re c ent in teraction with other cluster members and 
lMartin&Friedlil ( ll997t) did not find any particular nuclear activ- 
itjQ. NGC 4900 is part of a larger sample of galaxies observed with 
OASIS and represents a typical object for the presentation of our 
analysis technique; papers by Cantin et al. will follow to present 
the results for the other galaxies. 

Using the high resoluti on data from OASIS, the spectral syn- 
thesis code_J;AVALSB_JronijDio^^ Robert! ( |2006|) . and models 
from [Gonzalez Delgado et" "alj ( l2005h . we characterize the galaxy 
components in a consistent way and propose scenarios for the evo- 
lution of the galaxy central region. We have developed an iterative 
technique which has the advantage to separate a very young stellar 
population from an older one by taking into account the emission 
lines from the gas associated with the young population, the nebular 
and the stellar continuum from this young population, along with 
the stellar absorption lines and the continuum from an old popula- 
tion. This technique is different from the one used in other studies, 
which ofte n do not intend to study in detail all the stellar popu- 
lation s (e.g. lAUard et al.ll2006l : lGanda et al]|2007l : lstoklasova et al.] 
l2009h but where, for example, the contribution from a young pop- 
ulation is neglected even if emission lines are seen. 

This paper is organized as follow. The OASIS data and their 
reduction are presented in §2 along with the description of the 
method applied to measure spectral indicators and the synthesis 
models used to gather the information for the gas and the stellar 



1 In 1999. NGC 4900 hosted a supernova (SN 1999br) which was located 
on its outskirts, far from the region observed with OASIS. 
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Figure 1. OASIS acquisition image of NGC 4900. The image was obtained 
without a filter. The rectangle indicates the field covered by the lenslet array 
used in the OASIS spectroscopic mode. It samples a region 790 pc X 690 pc 
centered on the galaxy coordinates. 



populations. A demonstration on how the indicators are applied is 
done in §3 along with a description of the integrated spectrum. In 
§4 we present spatial maps of the continuum flux, emission lines, 
and extinction. A detailed analysis for all the spectra is done in §5 
where we first describe our iterative technique to separate in an ac- 
curate way the contribution from the different stellar populations. 
A discussion of the iterative technique and an evolutionary sce- 
nario for NGC 4900 are presented in §6. Finally, our conclusions 
are summarized in 57. 



2 OASIS DATA 

2.1 Observations and Data Reduction 

The data for NGC 4900 has been collected with the instrument OA- 
SIS at the 3.6 m Canada-France-Hawaii Telescope on March 18 
and 19, 2001. OASIS was used in the Cassegrain f/8 mode with 
the EEV detector and without adaptive optics. OASIS is an integral 
field unit which is composed of a matrix of lenslets. Documen- 
tation for OASIS may be found on the web page of the William 
Herschel Telescope, where it is now located. The spatial resolu- 
tion of 0.4 1 "/lenslet was selected to allow a large field of view of 
I5"xl2". Fig. [T] shows the OASIS acquisition image along with 
the location of the lenslet array used in the spectroscopic mode. 
Two wavelength ranges have been covered using the grating MRl 
(4760 to 5558 A) and MR2 (6210 to 7008 A). These configurations 
offer an instrumental spectral resolution of 2.17 and 2.23 Apixen^, 
respectively. 

One exposure of 2400 s was obtained for each of the gratings. 
Corresponding sky exposures were limited to 900 s due to time 
constraints. A set of calibration files (bias, dark, and micro-pupils) 
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Figure 2. The OASIS integrated spectrum of NGC 4900. The MRl configuration is shown in (a) and the MR2 configuration in (b). The features near 6280, 
6345, 6495, and 6850 A are caused by sky fine subtractions and telluric absorptions. 



along with spectra of Neon and Mercury lamps and a spectropho- 
tometric standard were collected to perform the calibrations. 

The data reduction was mainly done with the XOASIS pro- 
gram (version 6.0) developed by the TIGER team at the CRAL- 
Observatoire de Lyon. With XOASIS, we performed the usual oper- 
ations, i.e overscan removal, and dark current and bias sustraction. 
We also corrected for weak fringing in the detector using dome flat 
exposures. 

The next step in the data reduction, the mask extraction, was 
more crucial. The mask was created using XOASIS where a model 
of the telescope and instrument optical path was considered for 
each of the lenslets. The mask was applied to each image to pro- 
duce a datacube where the spectra were automatically (but roughly) 
wavelength calibrated. The refinement of the wavelength calibra- 
tion was done using the neon lamp image obtained in the same 
telescope position and instrument configuration as the science files 
to be calibrated. This gave a final wavelength resolution in MRl 
and MR2 of 2.3 ApixeP^. The datacubes were corrected for the 
instrument transmission using the flat field task. This was followed 
by the XOASIS routine to remove comic rays. The sky background 
was subtracted using an image of the sky in the case of the galaxy 
or using the image's surrounding sky information in the case of the 
standard star. 

The last steps done with the XOASIS software consisted in 
the flux calibration and the merging of the datacubes. We used the 
XOASIS flux calibration task. Considering the noise level and the 
uncertainty in the reduction process, we estimate the flux level to be 
accurate to within 10 per cent. Finally, we merged each calibrated 
datacube with itself. While this operation did not change the signal- 
to-noise ratio or any characteristics of the spectra, it simplified the 
analysis by reordering the lenslet identification number from a con- 
centric pattern to a rectangular system starting with lenslet #1 in the 
lower left South-East corner, and moving up North then West, col- 
umn by column. 

For the analysis, the datacubes were easily manipulated with 
IRAlQ. We corrected for the galactic reddening using the law 



Table 1. Line parameters 



Line 


Center 


Blue continuum 


Red continuum 




(A) 


band (A) 


band (A) 


H/3 


4861 


4841 


4881 


[Olll] 


4959 


4939 


5027 


[Olll] 


5007 


4987 


5027 


Mg2 


5159-5195 


5082-5140 


5210 - 5248 


Fez 


5270 


5250 


5290 


Fel 


5335 


5315 


5355 


[Ol] 


6300 


6280 


6320 


[Oi] 


6364 


6344 


6384 


[Nil] 


6548 


6523 


6603 


Ha 


6563 


6523 


6603 


[Nil] 


6583 


6523 


6603 


[Sii] 


6716 


6696 


6750 


[Sii] 


6730 


6696 


6750 



of ICardelli et"aD ( Il989l) and an extinction £fB-VJ=0.024 (NED) 
and then deredshifted the spectra using the galaxy velocity of 
960kms~^ (NED). To take into account the 0.8" seeing hmit at 
the time of the observation, we spatially smoothed over 4 adjoining 
lenslets for a final spatial resolution of ~ 53 pc (one pixel on the 
images will cover ~ 26 pc). Rejecting the low signal spectra on the 
edge of the image, this leaves us with an almost rectangular matrix 
of dimension 30x26 lenslets, or spectra, where the signal-to-noise 
ratio per pixel reaches 40 in MRl (measured between 5050 and 
5 150 A) and 50 in MR2 (between 6400 and 6500 A). 

2.2 Measurement of Spectral Lines and the Synthesis Codes 

The selected instrument configurations allow us to see a number of 
emission lines which will be used to characterize the stellar popu- 
lations, the extinction along with other gas properties, and nuclear 
activity, if any. These emission lines are: H/3, [O III] AA 4959, 5007, 



^ IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 



Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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[Ni]AA5200, [O I] AA 6300, 6364, Ha, [Nil] AA 6548, 6583, and 
[S II] A A 6716, 6730. In this wavelength range, we can also identify 
a number of stellar absorption lines which will be useful. These 
lines are Mg2 around 5177 A (which is composed of the magne- 
sium molecular band MgH and the magnesium absorption triplet 
Mg;,), Fe I A 5270, Fe I A 5335, and H/3 (when not contaminated by 
nebular emission). The list of the spectral lines measured is given 
in Table[T]along with the position of continuum bands on each side 
of the lines. An example of the spectral features observed may be 
seen in the integrated spectrum of Fig.|2] 

We estimated that for each invidivual spectrum, the galaxy ve- 
locity dispersion is small compared to the spectral resolution since 
each lenslet covers a relatively small spatial region. Therefore, we 
did not apply any correction to the spectra for this effect. Spectral 
line measurements were performed with an IRAF script to assure 
the uniformity of the result over the large number of spectra. The 
automatic measurements required that we first normalize the spec- 
trum continuum. In this routine, two continuum bands (see Table[T) 
on each side of a line must be specified. The lines were not mea- 
sured with the standard indicator definitions from the Lick system 
chiefly because our data has a higher spectral resolution than the 
one used in this system. A Gaussian was fitted to the line using the 
theoretical line center as a first reference for the Gaussian centroid. 
Most lines were fitted using a single Gaussian, but some (Ha and 
[Nil] AA 6548, 6583 or [Sll]AA 6716, 6730) needed a deblending 
routine. For each line, we extracted a center, an equivalent width 
(EQW), and a full width at half maximum. In the case of the Mg2 
line, a n index was calculated according to the definition of Worthe)! 
l ll994l) . A flux value was also calculated as the product of the Gaus- 
sian flux with the continuum value under the line as given by the 
non-normalized spectrum. Less than 30 per cent of the automatic 
measurements failed (giving peculiar values for the line parame- 
ters or extreme numbers compared to those obtained for the sur- 
rounding spectra) mainly due to a poorer signal-to-noise ratio in 
the continuum and a weaker line intensity. These have been redone 
manually. 

The H/3 emission was carefully measured at the base of the 
absorption profile when observed. No absorption component un- 
derneath the Ha emission was obvious and no correction was 
done. Prior to all calculation, special care was given to properly 
align the MRl and MR2 images. It was accomplished by super- 
posing the continuum peaks observed in the two configurations. 
A small shift, 1 pixel to the North and 2 pixels to the East for 
MRl relative to MR2, was applied. A total of 631 spectra are 
found to be overlapping in the two fields of view. The extinc- 
tion is then measured considering the ratio of the Ha and H/3 
e mission lines for a model as suming the recombination case B 
of lOsterbrock & Ferlandl ( |2006|) . The gas abundance 12+log(0/H) 
and the ionization parameter q are estimated using a m ethod 
similar to the one proposed by iKewlev & Dopit3 ( l2002h with 
the flux ratios: [N II] A 6583/[S II] A 6716+6730, [N II] A 6583/Ha 
[Nil] A6583/[OIII]A5007, and [O ill] A 4959+5007/H/?. The 
met allicity determinatio n depends strongly on the calibrators used 
(seelKewlev & EllisoiT| [2008). Other calibrators, such as those of 
IPettini & Pagell j2004h . available in the wavelength range con- 
sidered here, would give a metallicity estimate shifted down- 
ward by up ~0.2 dex compared to our result. In our work, 
we assume that solar met allicity is defined by 12+log(0/H)=8.7 
( iGrevesse&Sai^fTggih . AGN activity is looked for based on 
the line r atios: [N II] A 6583/H a an d [S II] A 6716-| - 6730/H a as pro- 
posed bv lDopita et"ZI j2006l) and [Baldwin et al.l ( Il98ll , hereafter 



BPT). The [Ol] A6300/Ha diagnostic is omitted since we cannot 
accurately measure the weak [O I] line in our spectra. 

The synthe sis ^code LAVALSB teobert et 2A\ l2003l : 

iDionne & RobertI 1 20061) is used to reproduce the emission 
lines EQW associated with very young star forming regions 
(^20Myr ). LAVALSB was de veloped in parallel with STAR- 
BUST99 ( Leithereret^ll999l) . It uses th e evolutionary tracks 
of the Genev a group, without rotation ( Schaller et al.^ 1 19921: 
iSchaerer et^ Il993jbl: Icharbon nel et al, 1993 ; Mevnet et al.l 
Il994h . at foiu^ metallicitie^ Z=0.002, 0.005, 0.02, and 0.04. For 
the young stellar populations we assume a metallicity identical 
to the one found for the gas. For this paper, the IMF used is a 
power law with a Salpeter slope a=2.35 and cutoff masses of 1 
and 100 M0. An instantaneous burst of star formation (where all 
the stars are formed at the same initial time) is considered. Nebular 
continuum flux, from the free-free and bound-free emission in 
hydrogen and helium, is also taken into account. The equivalent 
width of the hydrogen emission lines is calculated directly from the 
number of ionizing H I photons. These values are to be compared 
to the Ha and H/3 equivalent widths observed in order to obtain the 
age of the young stellar population. The flux level from the young 
population, including its nebular contribution to the continuum, is 
to be used to estimate the mass of the young population. 

Measuremen ts of the o bserved absorption lines are compared 
with models from lGonzalez Delg ad o et al . (2005) to study the stel- 
lar populations directly, loonzalez Delgado e t al. give synthesized 
spectra for metallicities of Z=0.004, 0.008, 0.019, and 0.030. These 
models allow the study of stellar populations with an age between 
4 Myr and 17 Gyr (or 60 Myr and 17 Gyr for the highest metallicity 
1.5 Z0). These models have been computed for a single burst pop- 
ulation, i.e. an instantaneous burst of star formation with no chem- 
ical evolution in the population beside normal stellar evolution. It 
considers an IMF with a Salpeter slope and masses between 0.1 
and 120 M0 . This smaller lower mass cut-off, compared to the one 
for young stellar populations, is required to properly replicate the 
oldest (>1 Gyr) population optical proprieties (Zhang et al.ll2002h . 
For populations younger than 20 Myr, the contribution to the op- 
tical flux from stars smaller than 1 Mq is negligible compared to 
the flux ofB and A type stars. The spectral resolution of the mod- 
els of Gonzalez Delgado et al. is high (0.3 A pixeP^) since they 
used the theo retical spectra for individual stars from the library of 
.Martins etal] ( l2005h . Therefore we degraded these models to match 
our observed spectra. Then, we measured the spectral indices over 
the synthetic spectra using the same method as for the observational 
data. From the absorption lines Mg2, Fe I A 5270, Fe I A 5335, and 
H/5 we obtain simultaneously the metallicity and age of a stellar 
population. The absolute flux from the population found is also 
used to est imate its contributio n to the toteil mass. O ther IMFs like 
the one of Kroupa et al.l ( Il993h and'Chabriejl ( l2003h would reduce 
the number of low mass stars in the models. This would slightly 
change the shape of the countinuum flux for the oldest populations 
and consequently modify the stellar mass of these populations. On 
the other hand, the age and metallicity values would probably re- 
main unchanged with these IMFs since these parameters are deter- 
mined through spectral indices that are mostly influenced by more 
massive stars. Testing other IMFs is beyond the scope of this paper. 



Here we adopt for the solar metallicity: Z0=O.O2 jAsplundetd]|2009h . 
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Table 2. Line measurements in the integrated spectrum 



Indicator 



Value 



EQW(Ha) 
EQW(H/3) 
EQW(Fe I A5270) 
EQW(Fe I A5335) 
l(Mg2 A5177) 
Log([N II] A6583/HO) 
Log([0 III] A5007/H/3) 
Log([S II] A6716+6730/Ha) 
Log([0 1] A6303/Ha) 



-38.5 ±0.1 A 
-8.2 ±0.2 A 
+0.9 ±0.1 A 
+0.4 ±0.2 A 
+0.065 ± 0.005 mag 
-0.46 ±0.01 
-0.40 ±0.01 
-0.60 ±0.02 
-1.5 ±0.2 




Figure 3. Equivalent widths of Ha and H/3 based on the synthesis code 
LAVALSB. The values for Ha (black curve) and H/9 (red curve) have been 
calculated for an instantaneous burst, at solar metallicity, and using a stan- 
dard IMF. The equivalent widths observed for the integrated spectrum are 
indicated as horizontal lines (their uncertainties are also considered, but are 
too small to be noticed on the plot). Both emission lines are consistent with 
a stellar population of ~ 7 Myr. 



3 INTEGRATED SPECTRUIM 

In order to guide us with the detailed analysis, to compare our data 
with previous studies from the literature, and to show how we use 
the stellar indicators and the models, we first present the integrated 
spectrum and its analysis. 

The integrated spectrum was created by adding together all 
the OASIS individual spectra. This spectrum covers a region with 
an average diameter of ^720 pc centered on the galaxy coordi- 
nates. The Balmer decrement reveals a global internal extinction 
£f £— Vj=0.25±0. 02 which fits very well with the published result 
of lHoetal.l ( fT997l) . The reddening corrected spectrum is shown in 
Fig. |2] The strong and narrow emission lines clearly indicate the 
presence of hot stars. The ]VIg2 and Fel absorption features are 
a more direct evidence for a stellar population, although without 
modeling, it is not easy to say if they belong to the same young 
population responsible for the emission lines or to an older under- 
lying population. 

Line equivalent widths and the spectral indice measured in the 
integrated spectrum are given in Table |2] In this analysis, all line 
measurements are used to obtain the stellar population character- 
istics independently of a possible contamination by other popula- 
tions. We notice that the nebular emission line ratios (see Table[2} 
do not indicate any peculiar activity for the integrated spectrum. 




Figure 4. The evolution of absorption lines with t he age and metallicity of 
the stellar population extracted from the models of lGonzalez Delgado et alj 
12005). The two top panels show the equivalent width of the iron lines and 
the bottom panel presents the Mg2 index. Different metallicities have been 
used as indicated in the top panel. The values observed for the integrated 
spectrum are indicated as horizontal lines (uncertainties are represented us- 
ing dash lines). 



These ratios are also in good agreement with those of IHo et aP 
l ll997l) . 

An oxygen abundance 12-flog(O/H)=8.70±0.05, i.e. Z0, is 
found for the gas and it is assumed to correspond to the abundance 
of the young stellar population responsible for the emission lines. 
The equivalent widths of the Ha and H/? emission lines are then 
compared with models from LAVALSB at this metallicity, as shown 
in Fig. [3] This comparison gives an age of 7.0±0.2 JVIyr for a single 
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Figure 5. Continuum and nebular emission maps. The continuum fluxes (10~^^ erg s~^ cm~^ A~^) from tlie MRl and MR2 configurations are shown in 
(a) and (b), respectively. The H/3 and Ha equivalent widths (A) are in (c) and (d), respectively. An equivalent width is given for the dominant component 
only, either the emission (negative value) or the absorption (positive value) component. Regions C 1 , C2, and C3 are delimited by contour lines extracted on 
the average continuum map. The region LI is defined by a contour line on the H/3 emission map. Only the area in which spectra are available for the two 
configurations is shown. On all maps North is up and East is to the left. 



young ionizing stellar population. The uncertainty is based on the 
equivalent width measurements using only a burst model at solar 
metallicity. For example, a small decrease in the model abundance 
to 0.4 Zq would only slightly increase the age of the young popu- 
lation and its uncertainty. 

Fig-Spresents the comparison of the observed measurements 
for the iron and magnesium absorption lines with the model values. 
These indicators should be equally used to get an age and a metal- 
licity. First, if we ignore the larger uncertainty for the Fe I A5335 
line (middle panel in Fig.|4l(, the three lines suggest a young pop- 
ulation with an age close to 7-8 Myr and a metallicity near solar 
(it is difficult to give a better estimate of the metallicity at these 
young ages as models with a metallicity higher that solar are miss- 
ing). This is close to the population found based on the emission 
lines. On the other hand, if a lower weight is given to the uncertain 
Fel A5335 line a much older population, between 120 to 630 Myr 
with a lower metallicity of 0.2-0.4 Zq, seems more likely. We did 
not use the absorption component of H/3 here because it is strongly 
contaminated by the nebular emission. In conclusion, the integrated 



spectrum clearly indicates that the galaxy central region is domi- 
nated by a metal rich 7 Myr population, but the characteristics of 
the population responsible for the absorption lines remain uncer- 
tain. 



4 MAPS OF CONTINUUM FLUXES, LINE EQUIVALENT 
WIDTHS, AND EXTINCTION 

In this section we take into account the spatial information to lo- 
calize the stellar populations and estimate the extinction in each 
spectrum. In the next sections, we will present a rigorous analysis 
of the stellar population parameters considering their relative con- 
tribution to the spectrum. 

Fig. |5(a)] and [5(b)l show maps of the continuum fluxes in the 
MRI and MR2 configurations, respectively. The flux was calcu- 
lated over the whole spectral range in each configuration, exclud- 
ing all features which would deviate from the average by 3 sigma 
or more. The maps display three peaks in the continuum flux. Con- 



tour lines around these peaks define the regions CI, C2, and C3. 
These are oriented from South-East to North- West along the large 
scale bar of the galaxy. The brighter peaks correspond to the two 
bright regions seen on the OASIS acquisition image (Fig. [T](. The 
region C2 corresponds to the location of the galaxy center. The im- 
portance of the flux in the region C3 is difficult to establish while it 
seems only partially included in the OASIS field of view. 

Fig. |5(c)| and |5(d)| are maps of the equivalent widths of the H/3 
and Ha lines. In these maps, the equivalent width of the emission 
component or the equivalent width of the absorption component is 
shown, depending on the strongest component in each spectrum. 
For example, while Ha is always in emission (a negative value for 
the equivalent width), an absorption component (positive value) is 
stronger in H/3 inside region C2. It is stricking to see the morpho- 
logical differences between the continuum and the emission line 
maps. Strong emission lines are seen near CI and C3. The emission 
line region LI is much more extended than CI, while its maximum 
is located ~80pc to the North of the CI continuum peak. Another 
emission region is also seen in the image South-East comer, right at 
the edge the field of view. The continuum region C2 barely displays 
any emission line; at this location, an H/3 absorption component is 
clearly found underneath a weak emission while Ha shows only a 
weak emission (see for example the spectrum in Fig.[9l(. Maps cre- 
ated for the other emission lines observed ([O III], [N II], and [S II]) 
are not presented here because they display the same morphology 
as observed for Ha. The emission lines are most likely the conse- 
quence of the interstellar gas ionized by the ultraviolet flux from 
OB stars present in young stellar populations. The continuum and 
emission line maps (not like the integrated spectrum) offer the best 
evidence for the presence and the superposition of both old and 
young stellar populations in the galaxy center. 

From the ratio of the emission lines Ha and H/3 (with a care- 
ful measurement of the line emission component only) important 
spatial variations of the extinction are found, as shown in Fig. |6] 
It indicates E(B—V) values between 0.02 and 1.3, with an average 
extinction of 0.47±0.08. Larger values are found in isolated pix- 
els, where the signal is low or the H/J emission line is very weak. 
For a few spectra (mainly in the North-East corner of the image), a 
clear H/3 emission component could not be observed over the noise 
and we adopted an extinction value computed from the average of 
four adjoining spectra, assuming that the extinction changes slowly 
from one spectrum (26 pc) to the next. The extinction map shows 
a peculiar morphology; a clear dust lane is seen superimposed to 
the continuum region C2, and oriented perpendicular to the galaxy 
large scale bar. This heavily dust-obscured zone, identified XI on 
the map, displays an average extinction of 0.62±0.06. Over the 
South-East boundary of this zone, near CI and LI, the region X2 
identifies a much lower extinction zone with an average E(B—V) 
of O.I9±0.01. The region X2 is surrounded by a background with 
an average extinction of 0.31±0.01. Over the North-West bound- 
ary of the dust lane, near C3, the extinction has an average value of 
0.39±0.01. 



5 CHARACTERISTICS OF THE STELLAR 
POPULATIONS 

5.1 The Iteration Technique 

For a consistent analysis of the stellar populations, we aim at disen- 
tangling the fluxes associated to each population before determin- 
ing their characteristics. An iterative technique is proposed in order 
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Figure 6. Extinction map. Only the area which overlap between the two 
configurations is shown (i.e. the H/3 line is part of the MRl configuration 
while Ha is in MR2). Areas of high (XI) and low (X2) extinction are iden- 
tified. Regions CI, C2, C3, and LI are as in Fig. [3] 

to isolate the different stellar populations present in a spectrum. The 
emission lines of H/3 and Ha, when present, are used to character- 
ize a young population (<20Myr) ionizing the nebular gaz, while 
absorption lines are used independently to describe directly a stel- 
lar population. The absorption lines may show the same age as the 
one described by the emission lines but, based on our experience 
with this galaxy and the others in our sample, the absorption lines 
mostly describe an older generation of stars than the one found with 
the emission lines. So for simplicity's sake, we say through the pa- 
per that we find a young population based on the emission lines 
and an old population based on the absorption lines. In any case, if 
only a young population is present the two sets of indicators will 
provide comparable characteristics for the population. Our itera- 
tive technique is summarized in Fig.|7l We use a program based on 
Bayesian statistics to compare line measurements between obser- 
vations and models. This program looks at a given set of indicators 
in order to identify a single age-metallicity pair solution or a single 
mass solution, with its uncertainty, as described in the Appendix. 
Our iterative technique begins with a preliminary analysis, named 
"iteration 0", which gives a starting point for the characteristics of 
the stellar populations. Iteration looks at the emission lines and 
the absorption lines without taking into account the contribution of 
the other population to the spectral lines studied (as done for the 
integrated spectrum). 

Iteration I consists in subtracting a theoretical model for 
the old population in order to study th e emission lines alone. 
This m odel is adapted from the models of lOonzalez Delgado et al.l 
( l2005l) . as described in §2.2, and is selected using the age and 
metallicity found after iteration for the old population. The con- 
tinuum flux level of this theoretical model (which is related to the 
mass of the old stellar population on the line of sight for a specific 
spectrum) is set such that the continuum flux left after its substrac- 
tion from the observed spectrum is as close as possible to zero with- 
out being negative. The spectrum resulting from the subtraction is 
then reprocessed as if from a single young burst alone. Bayesian 
statistics are used to compare line measurements between obser- 
vations and LAVALSB models. To complete iteration I, the model 
which represents the old stellar population and the best model for 
the young population are summed. A nebular continuum calculated 
using LAVALSB based on the ionizing flux of the young population 
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Figure 7. Schematic representation of the iterative technique. 



is also included. The sum of these models are then scaled so that 
when compared with the observed spectrum (using again Bayesian 
statistics), it reveals a new mass for both, the old and the young, 
stellar populations. 

Iteration 2 consists of subtracting the theoretical model for the 
young population (including the nebular continuum) found in the 
previous iteration from the observed spectrum in order to study the 
absorption lines. Prior to the subtraction, synthetic emission line 
profiles of Hof and H/3 only (no synthetic lines are used for [O III], 
[Nil], or [S II]) are added to the spectral energy distribution of the 
young population spectrum. It is assumed that the Hq line profile is 
the same as the one observed (we consider that the absorption com- 
ponent underneath Ha is negligible), while the H/3 line is a copy 
of the Hq profile scaled with respect to the theoretical 2.86 line in- 
tensity ratio used for the extinction correction (see ^2.2\ . Once this 
theoretical spectrum is subtracted, revised values for the age and 
metallicity of the old stellar population are obtained. At this step, 
the H/3 line might still present some emission and cannot be used 
for determining the characteristics of the old population. As the it- 
eration proceeds, the H/3 profile should present a pure absorption. 
The masses of both stellar populations are recalculated by compar- 
ing the combined model for the two populations with the observed 
spectrum. 

The next steps consist in repeating iterations 1 and 2 until there 
is no significant change, from the previous step, in the age, metal- 
licity, and mass found for both populations, or until the difference 
for the stellar indicator measurements are of the order of their un- 
certainties. In general for most spectra, we need to repeat iterations 
1 and 2 three times. In the following subsections, we present the 
results from iteration and, in more details, the characteristics ob- 
tained after the final iteration. 



5.2 Iteration 

We corrected the individual spectra with their respective extinc- 
tion value (as given in Fig. [6} using the reddening law of Cardelli, 
Clayton, & Mathis (1989) within IRAF. As done for the integrated 
spectrum, we estimated the age and metallicity of a young popula- 
tion, based on the emission lines, and an old population, based on 
the absorptions. This analysis is done again without taking into ac- 
count the contribution to the continuum from one population when 



studying the other. It will be used as a starting point for the iteration 
process. 

At this step in the iteration process, most of the lenses in the 
observed region show emission line ratio typical of H II regions 
(only a few spectra, notably in C2 and on the periphery, display 
emission line ratio values as found in AGNs). The emission lines 
reveal a metallicity for the gas, measured in term of 12-l-log(0/H), 
which is slightly above the solar value along the C1-C2-C3 axis, i.e. 
along the galaxy large scale bar. The gas metallicity remains closer 
to the solar value away from this axis (but in these regions it is also 
more difficult to make good estimates because of the lower signal 
in [O III] A 4959). Adopting the gas metallicity for the stars respon- 
sible for the ionizing flux, we find that the youngest regions are 
located in the South-East image comer (see Fig. |8(a)[ l. We obtain 
an age of about 5.7±0.3 Myr at the location of the stronger emis- 
sion peak in LI, an average age of 6.5±0.3 Myr in regions CI and 
C3, and 12±1 Myr in C2. The age given by the Mg2 and Fel lines 
(see Fig. lSMt , 

is on average over the whole map, 750±250 Myr. 
There is no clear spatial region with a different and specific age ex- 
cept for a few patches which may look younger (for example near 
CI). The most probable metallicity given by the absorption lines is 
Z=O.2Z0. 

5.3 Final Iteration 

Using the results from iteration 0, we follow the steps described for 
the iteration technique (Fig.|7}. The final results for the age, metal- 
licity, and mass of the young and old populations are presented here 
along with examples of the final model for selected spectra. As de- 
scribed at the end of this section, the iteration technique also helps 
us to reveal a low level of peculiar activity in the galaxy center. 

Fig. m shows the best model obtained after the final iteration 
for three spectra as examples. One of them is located in CI, where 
the young population is dominating the flux, another one is located 
near the galaxy center, in C2, where weak emission lines are seen, 
and the last spectrum was taken North of CI, in XI, where both a 
young and an old population are more or less equally contributing 
to the flux. This figure presents the individual contributions to the 
final model from the young and the old components, along with the 
residual between the observed spectrum and the final model. Since 
the best models are found using measurements in many lines and 
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(a) 




Figure 8. Age distributions. At the top are maps obtained after the iteration 0. 
in Myr, based on the emission lines. To the right, (b) and (d), the age is given, 

in the continuum, a perfect superposition is not easy to achieve. 
The overall residual is generally close to zero (and will never be 
equal to zero in the emission lines of [O III], [N II], and [S II] since 
they are not taken into account in the young population models) 
and indicate a good fit of the final models. The spectrum in CI 
(top panels in Fig.|9j is best reproduced by a young population of 
6.8±0.2 Myr with a metallicity of 2 Zq (based on an oxygen abun- 
dance of 9.1±0.3) and a mass of 3.9±0.4 x 10^ M0, superposed 
with an averaged old stellar population of 850±150Myr with a 
metallicity of O.2Z0 and a mass of 3±1 xlO^M©. In this case, 
the flux from the young stellar population is so important that it 
is not possible to constrain the characteristics of the old compo- 
nent. Therefore, an averaged old population spectrum was simply 
adopted. In C2 (middle panels in Fig.|9]l, the best model was found 
by co-adding a young population of 7.8±0.2Myr at 2Zq (based 
on an oxygen abundance of 9.0±0.2) and a small stellar mass of 
3.5±0.4 X 10'' M0, with an old stellar population of 180±50Myr 
at 0.2 Zq and a mass of 7.8±0.5 X 10^ Mq . The best model for the 
spectrum in XI (bottom panels in Fig.|9]l was obtained by adding 
a young population of 7.5±0.3Myr at 2Zq (based on an oxygen 
abundance of 8.9±0.4) and a very low mass of 6.0±0.5 x 10* Mq, 
with an old stellar population of 3±1 Gyr at 0.47© and a mass of 
4.0±0.3 xlO^ Mq. 



(b) 




At the bottom, after the tinal interation. To the left, (a) and (c), the age is shown, 
ing a logarithm scale, by the absorption lines. Contour lines are as in Fig.|5] 

5.3.1 Age and Metallicity 

Fig. [8] also shows two maps of the distribution of the age for the 
stellar populations after the final iteration, one based on the emis- 
sion lines (Fig. |8(c)[ l which is limited to younger populations and 
the other based on the absorption lines (Fig. [8(d)}. As found after it- 
eration 0, younger populations are located near regions CI and C3. 
The age at the location of the emission peak in LI is 5.5±0.1 Myr 
and the average age is 6.0±0.1 Myr in CI and 6.4±0.1 Myr in C3. 
In C2, a few spectra display ratios of the emission lines that indi- 
cate some peculiar activity in the nucleus (see § 5.3.4). Excluding 
these spectra, an age of 7.9±0.3Myr is found in C2 based on the 
emission lines. Surrounding C2, along the dust lane, the emission 
lines give and average value for the age of 8.2±0.4Myr. Note that 
in the upper North-East comer of this map, the age is often given 
by the Ha line only as an emission in H/3 is not always easy to 
measure. 

The absorption lines reveal a more or less homogeneous old 
population of 850±150Myr, on average, distributed all over the 
central galaxy region observed (Fig. |8(d)[ l. Within the region C2, 
the average age is somewhat older, 940± 120 Myr, but no other spe- 
cific pattern or gradient is seen for the old population age. It is dif- 
ficult to study the absorption lines near CI because the contribution 
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Figure 9. Examples of best models for three different regions. Top panels: spectrum # 322, in CI, dominated by strong emission lines. Middle panels: spectrum 
# 496, in C2, dominated by an old stellar populations. Bottom panels: spectrum # 330, in the high extinction zone XI at the North of CI, where the flux from 
the two stellar populations is nearly equals. The black line represents the observed spectrum in the MRl configuration (left) and the MR2 configuration (right). 
The contribution of a young stellar population characterized using the emission line is shown by the blue curve while the old stellar population studied using 
the absorption lines is represented by the red curve. The best model obtained by combining the contribution from the two stellar populations is shown by the 
red dotted curve. The blue dotted curve is the residual obtained after subtraction of the best model from the observed spectrum. Characteristics of the stellar 
populations used for the models are indicated on the right panels. 



to the continuum from the young population is very strong. After its 
subtraction, the signal-to-noise ratio in the absorption lines is very 
low and models give unreliable values (often the lower age limit 
from the models with an uncertainty equals to zero or to an ex- 
tremely large value). Therefore, in the map based on the absorption 



lines, we have replaced unreliable values (due directly or indirectly 
to a low signal-to-noise ratio) by white pixels. 

The final iteration reveals values for the age of the stellar pop- 
ulations which are significantly different from those obtained after 
iteration 0. A comparison of the age maps in Fig. [8] clearly shows 
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Figure 10. Metallicity distributions after tlie final iteration. In (a), tlie gas 
metallicity 12+log(0/H) is sliown. In (b), the stellar metallicity Z is given. 
Isolated pixels with a lower metalhcity value in (a) display very week emis- 
sion Hnes. White pixels in (b) indicate poor line measurements because of a 
low signal-to-noise ratio for the absorption lines studied. Contour lines are 
as in Fig. [5] 

this. Numbers for the young population age are slightly but system- 
atically smaller after the final iteration mainly because the stellar 
absorption component underneath the emission lines is subtracted 
in the iteration process. The old population age is, on the other 
hand, increased after the final iteration. This increase is brought 
forth mostly by two effects, first because the subtraction of the 
young population continuum flux has a tendency to increase the 
equivalent width of the absorption lines and second, because the 
addition of the Hfi absorption feature to the age determination pro- 
cess results, in this case, in a larger probability for an older age 
value. 

There is not a significant difference between the metallicity 
calculated after the final iteration and iteration 0. The gas and stel- 
lar metallicity maps are almost featureless (Fig.llOt. After the final 
iteration, the average gas metallicity is supersolar at 12-l-log(0/H)= 
8.9±0.1 (i.e. 2Z0) while the stellar metallicity is lower with an 
average value of Z=O.2±O.2Z0. Inside the region CI, the gas 
metallicity increases to an average value 12-l-log(O/H)=9.1±0.1, 
i.e. 3 Zq . For the same region, the estimated stellar metallicity is 
unreliable because of the very strong flux from the young popula- 



tion which reduces the signal-to-noise ratio in the absorption lines. 
In C2, if we assume that the peculiar activity seen in a few spectra 
has a negligible effect on the continuum and absorption lines, we 
find on average a larger stellar metallicity, Z~ 1 .0±0.2 Z© . A larger 
value for the stellar metallicity in the galaxy center is not surpris- 
ing. The difference between the stellar and gas metallicity is, on the 
other hand, quite important. 

To get a better idea of how the various populations found are 
distributed, we plot in Fig. [TT] two histograms showing the num- 
ber of spectra as a function of the age. The younger populations, 
found using the emission lines, are distributed in a narrow spike, 
centered at about 7 Myr (Fig. |l l(a)^ . An asymmetry indicated by a 
second spike, much weaker, near 6 Myr, corresponding to the av- 
erage population in C 1 . The histogram based on the age extracted 
from the absorption lines shows an old and long episode of star for- 
mation in the galaxy central region (Fig. |1 l(b)| (. It started about 
3 Gyr ago, was very active around 1 Gyr, and lasted until about 
100 Myr. There is also an indication of a star-forming activity be- 
tween 6 and 30 Myr and around 10 Gyr in this histogram, but it is 
associated with sporadic pixels on the age map which are often at 
the limit of our rejection criteria due to a lower signal in the ab- 
sorption lines. 



5.3.2 Mass 

Fig.[T2]shows the mass distributions for the young and the old stel- 
lar populations, individually, as well as for the total stellar mass 
derived in each spectrum. The mass of the stellar populations is 
obtained during the iteration process by comparing the continuum 
level of the observed spectrum with the best model selected based 
on the line indicators after each time step. The mass estimated 
here are current values. The current mass of a stellar population 
is smaller than the initial mass as a consequence of stellar evolu- 
tion; the most massive stars evolve and die more quickly and do not 
contribute to the flux and mass observed now. We estimate, based 
on models for an instantaneous burst with a standard Salpeter IMF, 
that the current mass of a 1 Gyr population would be ~26 per cent 
smaller then the initial mass, and for a 10 Gyr population it would 
be ~40 per cent smaller. We find the highest mass concentration 
of the young populations along the axis C1-C2-C3 and underneath 
the dust lane (Fig. |12(a)} . In CI, where the young population is very 
important, we have considered the average 850 Myr old population. 
If this population was younger, 10 Myr for example, then the mass 
of the young population in CI would not change significantly, i.e. 
by less than a factor of 10. Typically in a spectrum in CI and C2, 
there is about 5±1 x 10'' Mq, corresponding to an average density 
of about 50 M0 pc~^ and 40 M© pc~^, respectively. Outside these 
regions there is an average mass of 8± 1x10* Mq per spectrum and 
an average density of 30 Mq pc^^ (except to the North-East where 
it is smaller by a factor ~ 10). Summing for all pixels, the total 
mass in new stars is 1.14±0.05 x 10* Mq. 

The old stellar population is more massive at the center than 
outside (Fig. |12(b)| , with values from 2±1 x 10^ M0 per spectrum 
in the region C2 (~ 20 Mq pc~^) to 8±3 x 10"* Mq per spectrum 
in the surrounding. The mass of the old population near CI is often 
not shown on the map. If we compute the mass in CI using an old 
population with the average age of 850 Myr, we find the same value 
as in the surrounding region. (If we consider instead a younger age, 
10 Myr for exemple, the mass found would be reduced by a factor 
of ~10). Summing for all pixels available, we find that the old stel- 
lar population represents 2.07±0.06 x 10* Mq, about 60 per cent of 
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Figure 11. Histograms of the number of spectra as a function of tlie age 
after tlie final iteration. In (a) the age is shown based on the emission lines. 
In (b) the age is given by the absorption lines (as for Fig.[8). 



the total stellar mass. It clearly shows the importance of the most 
recent star forming episode that occurred about 7 Myr ago. 

Fig. |12(c)| shows the map for the total stellar mass obtained 
after adding together the two stellar populations characterized with 
the iterative technique. The mass is strongly concentrated in the 
galaxy center in C2. The region CI, where we are considering 
the average underlying older population, is also important. This 
is similar to the continuum maps (see Fig. |5(a)| and |5(b)| (. We 
also distinguish the galaxy large scale bar in the total stellar mass 
map. Within the section of the bar observed (supposing a width 
of 260 pc), the average density is 40±5 pc~ and it includes 
2.2±0. 1x10* Mq , i.e. almost 2/3 of the mass present in the galaxy 
central area. The region C2 contains about half of the mass of the 
bar observed. 



5.3.3 Peculiar Activity 

Fig. [m presents maps of the emission line ratios 
log([Nll] A6583/HQ) and log([S II] A 6716+6730/Ha). In 
general, both maps show typical ratios for Hll regions, i.e. 
where [Nil] and [Sll] are weak compared to Ha. The dust lane 
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Figure 12. Mass distributions after the final iteration. The maps show the 
current stellar population mass, in log(mass [M©]). In (a), the mass of 
the young population based on the emission lines. In (b), the mass of the 
old population based on the absorption lines. The total population mass is 
shown in (c). White pixels in (b) indicate missing measurements because 
of a low signal-to-noise ratio for the absorption lines studied. Contour lines 
are as in Fig.|5]and|6] 



Stellar Populations in NGC 4900 13 




I 



-0.3 



-0.4 



-0.6 



-0.8 



-0.4 



I 



-0.81 



Figure 13. Non-thermal indicator maps. In (a), log{[N II] A 6583/Hq:) and 
in (b), log([S II] A6716-l-6730/Ha). Contour lines are as in Fig.|5]and[6] 



(region XI) displays ratios which are closer to zero, but this 
is where the emission lines are weaker and, particularly in the 
North-East image corner, the signal is also lower. On the other 
hand, the region C2 has good signal and it shows a ratio of 
log([Nll] A6583/Ha)>-0.2 which is hig h for an H II reg i on ari d 
more typical for an AGN according to IStasinska et al. l l2006h . 
However the value of log([S II] A 6716+6730/Ha) in C2 remains 
typical of an H II region. 

The BPT diagrams, generally used to identify the type of the 
activity present in a galaxy, confirms the absence of a strong non- 
thermal activity. In Fig.[T4]the red solid line the top panels indicate 
the li mit between H I I regio ns and Seyfert/LINER activities accord- 
ing to lKewlev et al] ( 1200 ll) . The dotted line in panel (a) defines the 
limit w here H II regions are observed according to Kauffmann et al.l 
( l2003h . The blue line in panel (b) separates Seyferts and LIN- 
ERs dKewlev et al]|2006h . No spectrum is lying above the AGN, 
Seyfert, or LINER limit. Nevertheless, we see a significant num- 
ber of spectra falling into the composite/transition region, i.e. be- 
tween the solid and dotted lines of the log([0 III] A 5007/H/3) vs 
log([N II] A 6583/Hq) diagram. Fig. |14(c3] shows the location of 
these spectra on the image. Most of them are near the galaxy center 
in C2. A few isolated spectra are located along the dust lane, where 
lower signal have been seen in the North-East corner. We therefore 




Figure 14. Activity indicators. The BPT diagrams of (a) 
log([0 III] A 5007/H/3) vs log([N II] A 6583/Ho) and (b) 
log([0 III] A 5007/H/3) vs log{[S II] A 6716+6730/Ha) are shown. 
Lenslet # are given on those plot. The map in (c) indicates the location of 
pixels where an emission associated with a composite/transition region is 
observed according to the diagram in (a). Contour lines are as in Fig. \5\ 
and|6] 
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conclude that NGC 4900 hosts a composite/transition nucleus, i.e. 
a nucleus where both star formation an d an AGN are present (as 
defined bv lWang & Weill2008l : lHdl2008h . 



6 DISCUSSION 

6.1 The Iteration Technique 

The iteration technique has allowed us to refine the evaluation of 
the stellar population age, metallicity, and mass within a two pop- 
ulation scenario. The age and mass of the old stellar population, 
based on the absorption lines, are the parameters that have im- 
proved the most. With the iterative process, the nebular emission 
in H/? can be removed, allowing us to use the H/3 absorption com- 
ponent as an additional diagnostic and to reduce considerably the 
error bar of the old population age. By itself, the age value has lit- 
tle influence on the final mass of the old population. On the other 
hand, the presence of an important young stellar population affects 
the continuum level which also changes significantly the determi- 
nation of the old population mass and age. The young stars con- 
tinuum emission, and its associated nebular continuum, dilute the 
absorption lines of the old stellar population. In the region C2, 
for example, we find a typical age for the old population around 
750±250 Myr before the iterative process, while after, the age in- 
creases to 940±120Myn The change is less spectacular for the 
young population, since the age and its uncertainty are, in this case, 
smaller (typically by 0. 1-0.2 Myr). Nevertheless the young popula- 
tion age is also improved as a result of the better treatment of the 
underlying absorption in the nebular emission lines and the contin- 
uum level. This implies that emission line intensities for H/3 and 
Hq are stronger than before and, in consequence, the age estimated 
from these indicators is younger. 

The stellar metallicity has not changed significantly after the 
iterative process. It is strongly dependent on the Fe I absorption 
lines, but mainly, it depends on statistics based on only a few large 
metallicity bins in the models. The oxygen abundance, however, 
is estimated from line intensity ratios, involving the H/3 and Ha 
lines which are better constrained as the iteration process converge. 
The influence is keenly seen in the [O III] A 4959+5007/H/3 ratio, 
although this is only a small part of the calculation, slightly in- 
creasing the oxygen abundance. 

It is also important to realize that with the use of an in- 
tegral field unit like OASIS we are able to resolve smaller star 
forming regions in nearby galaxies. The stellar population in one 
lenslet may become so small in mass that we could be losing 
the statistical aspect of t he ini tial mass function. According to 
ICervifio & Valls-Gabau j ( |2003|) a mass below lO'^Mo, particu- 
larly for a young population, is becoming critical. On the other 
hand, we think that the statistical study of characteristics found over 
many lenslets is compensating for this. 

The iterative technique allows the identification of a weak pe- 
culiar activity in the central region of the galaxy. Before the itera- 
tive process, the underestimation of the H/3 absorption flux tended 
to assign a non-thermal activity to many lenslets, even in some 
located far away from the galaxy center. With the iterative pro- 
cess we have refined the line ratios to a point where most values 
are typical of H II regions while only a few spectra, mostly those 
concentrated in the galaxy central region, belong to the compos- 
ite/transition region of the BPT diagram. In the case of NGC 4900, 
this non-thermal activity is very weak and did not cause difficulty 
for the analysis. A stronger non-thermal component is expected to 



affect the continuum level, the absorption lines, and also the extinc- 
tion estimation. The presence of a strong non-thermal component 
will introduce an uncertainty for the population parameters with an 
iterative technique that only takes into account the superposition of 
stellar populations. 

Two aspects that have been neglected in our iterative tech- 
nique, and deserve a better treatment in the future when possible, 
concern the young stellar population contribution to the H/3, Mg2, 
and Fe I lines and the determination of the extinction. If the metal- 
licity of the young stellar population is below solar and the old 
population is older than 1 Gyr, it is not really a problem to neglect 
the absorption lines created by the young population since the old 
population will completely dominate these structures. However, a 
young population with an age between 10-30 Myr and an abun- 
dance near or higher than solar may contribute to these lines, as 
can be seen in Fig. |4l and influence the characteristics found for 
both stellar populations. For the moment, the models we are using 
do not include the absorption lines in all young stellar populations 
and therefore we could not consider this effect in our study. The ex- 
tinction was estimated in a simplified way in our study. It was done 
in the early phase of the analysis, prior to the iteration process, us- 
ing the overall Ho? line (without any correction for an absorption 
component) and our best measurement for the emission component 
in H/? (at the very base of an absorption component when present). 
A specific extinction was calculated for each spectrum, allowing us 
to draw an extinction map, which is certainly a big advantage for 
our study. Another improvement would be to add a treatment of the 
extinction into the iterative technique and see its effect on the stel- 
lar population parameters. As the emission lines get better defined 
along the iteration process (i.e. as the underlying stellar absorption 
components are better known), it will modify the extinction, the 
slope of the continuum, and also influence the choice of the best 
model for the stellar populations. A preliminary test indicated, be- 
side a much longer iteration time, a modification of the mass of 
the stellar populations (but not more than by a factor of 2) with 
a negligible effect on the age, gas oxygen abundance, and stellar 
metallicity. 

6.2 An Evolutionary Scenario 

The information collected provides clues about the history of the 
central region of NGC 4900. With these clues we propose here a 
simple evolutionary scenario leading to the structures and popu- 
lations observed. But first, let us summarize and further analyze 
these clues. We find a few young bursts, with an age between 5.5 
and 6.5 Myr, located along the galaxy large scale bar. We see an ex- 
tended region of about 8 Myr superposed to a dust lane which is ori- 
ented perpendicular to the galaxy large scale bar. These young pop- 
ulations are superposed to a longer episode of star formation with 
an average age near 850 Myr. In the galaxy nucleus, the age of the 
old population reaches 940 Myr and there is a composite/transition 
region activity. Furthermore there is a clear distinction between 
the gas metallicity, 27.Q on average (and higher near one of the 
younger bursts), and the old stellar population sub-solar average 
metallicity. 

The stellar mass associated to the episodes of star formation 
shows an interesting distribution. Fig. [15] shows histograms of the 
initial mass for each populations as a function of age. As found for 
the current mass, nearly 40 per cent of the total mass converted into 
stars took place during the most recent bursts. This star formation 
was also very efficient, i.e. 83 per cent of the stellar mass in these 
bursts was created between 6 and 8 Myr Most of these young stars 
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Figure IS. Histograms of the stellar population mass as a function of age. 
The initial mass is represented here. In a), for the young populations based 
on the emission lines and in b), for the old populations based on the absorp- 
tion lines. Different colors are used to show the mass in each age bin as a 
function of the stellar population metallicity (as in Fig llO) . 



are showing a metallicity that is twice solar, while the remaining 
young stars have a higher metallicity (3 Zq ) and dominate the mass 
budget bellow 6 Myr. The old star-forming episode, responsible for 
about 60 per cent of the inital stellar mass, looks very different. We 
interpret this episode as a long series of multiple bursts which took 
place in the past in the central region of the galaxy. The amplitude 
of the star formation was variable with time, reaching a maximum 
about 1 Gyr ago. Furtheremore, the stellar metallicity was changing 
with time for these old stellar generations as well. At the 1 Gyr 
peak, an important part of the old populations had reached the solar 
metallicity while before, the metallicity was lower. It is interesting 
to see that among the more recent old generations (after the 1 Gyr 
peak), the metallicity is lower than solar. 

Using the mass histograms, we also built a plot of the total 
stellar mass created before a certain age. The first derivative of this 
mass may be interpreted as the star formation rate at a particular 
moment in the history of the galaxy. Looking at Fig.[T6l we find that 
the historical star formation rate was small (below 0.3 Mq yr~^) 




Figure 16. Historical star formation rate. In black, the total mass formed 
before an age, and in red, the star formation rate for this age. In a) for the 
whole age range studied here (from 4 Myr to 12 Gyr), b) a zoom for the 
period covering the last 2 Gyr, and c) the most recent period between 4 and 
50 Myr. 



until very recently, when the most recent star-forming episode 
shows a rate of almost 40 M© yr~^. This rate is much higher then 
what would be found using the Hq luminosity of the integrate d 
spectrum (0.21 yr~^ with the equation of lKennicutllll998h . 
Before the recent 6-8 Myr burst, we see 3 bumps in the historical 
star formation rate: around 150 Myr, near 500 Myr, and at 1 Gyr 
corresponding to the maximum in the mass distribution of the old 
populations. We must point out that this method probably underes- 
timate the historical star formation rate since the width of the age 
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bin increases with time. If we consider tliat the width of the most re- 
cent burst {^li Myr) is typical for an episode of star formation, the 
historical star formation rate for the older bumps would increase to 
7M0 yr~^, 16Mq yr~^, and 40 yr~^, respectively. This last 
number is comparable to the SFR of the most recent star forming 
event. 

The narrowness of the histograms for the young stellar popu- 
lation mass (Fig. |15(a)| l or the number of spectra (Fig. |l l(a)| ( as a 
function of age, obtained based on the emission lines, may suggests 
that the young bursts have been triggered by a single specific event 
(internal or external) in the recent past. This event had an influence, 
to some level, on the galaxy bar which then became more active 
in compressing the interstellar material and triggering or fuelling 
new star formation, as proposed by different simulations (see for 
example|Athanassoulal 1992l ; lNorman et al.ll996l : lRegan & TeubenI 



l2004 l Haan et al. .I l2009lr The fact that we see young bursts located 
along the galaxy bar supports the importance of the bar. Enriched 
interstellar material may be the result of past episodes of star for- 
mation seen all over the galaxy central region. However, the fact 
that we see a higher oxygen abundance for the gas near one of the 
recent burst (in the region LI) and the fact that we find a large 
amount of material in this region support a flow of enriched ma- 
terial toward the central region. We are not seeing star formation 
in a nuclear ring, but its radius could be over the limit of our field 
of view (which has an average di ameter of ~720 pc com pared to 
a typical ring diameter of 1 kpc; iButa & Crocked [l 9931) . On the 
other hand, the extinction map shows clearly a dust lane oriented 
perpendicular to the galaxy large scale bar, which may be inter- 
preted as a dusty nuclear bar. A nuclear bar is also often cr eated 
in large scale bar simulations tenglmaier & Shlosmanll2004h . Nu- 
clear bars and star form ing knots are signatures of secular evolution 
jKormendv & Kermicutt 2004) which is the mechanism used to ex- 
plain the formation of pseudobulge in late spirals. 

The velocity map drawn using the centroid of the stronger 
emission line of Ha is not particularly clear with respect to the 
idea of gas from the disk flowing toward the galaxy center, but it 
surely indicates important motions. This map (Fig.|17b shows vari- 
ations in velocity {v sin i, where the galaxy inclination is of the 
order of 20°) from —300 to 200kms~^ with some clear patterns: 
the central region has a velocity close to zero, the South-West edge 
of the image displays a redshifted velocity, while near and along 
the large scale bar we see blueshifted regions not necessarily su- 
perposed to the young star forming knots. One may see on the East 
side of the dust lane a redshifted velocity, while the West side dis- 
plays mostly a blueshifted velocity. But overall this map is not easy 
to interpret. Velocity maps for the other lines (e.g. H/3 in emission 
and absorption, and Mg2; not shown here) are significantly more 
noisy because of the line w eakness. No ro tation curves have been 
published for NGC 4900, but lWarmelsl |T%8) gives a maximum ve- 
locity of 1 92 km s ^ ^ , based on the H I observations, which describe 
the motion in a typical disk, far away from the galaxy center. 

The old star formation episode appears to be more or less uni- 
formly distributed over the central region with older and more metal 
rich objects closer to the nucleus. This event is distributed over a 
long period of time; it took place between 100 Myr and 3 Gyr ago. 
It may indicate again the past signature of slow secular evolution 
which contributed to an initial phase of the pseudobulge forma- 
tion. Was it purely internal? The amplitude of this old star forming 
episode, larger that the more recent ~ 7 Myr episode, may suggest 
that it has a different origin although the star formation rate are 
possibly comparable. Was it related to an int eraction with another 
member of the Virgo cluster (see for exemple: lElizondo et al .1 19991 : 
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Figure 17. Velocity map (ii sin « in km s ^ ) based on the Ha emission line 
center. Contour lines are as in Fig.|5]and|6] 



iBoselli & Gavazzil2006l) . or did it involved a larger reservoir of gas 
already available at the time in the disk of the galaxy? In a simple 
scenario, all episodes of star formation could just be the conse- 
quence of the galaxy large scale bar evolution. If process such as 
bar buckli ng (i.e vertical instabil ity that thickens and weakens the 
system bar lOebattista et aill2004h combi ned with the tri ggering of 
star formation from infall from the bar i' Woznial3l2007 *). could be 
repeated it could produce many star formation episode in the galaxy 
center reproducing the populations and structures observed here. To 
verify this hypothesis, it would be useful to apply our technique to 
more regions on the galaxy bar and disk at different radius, to look, 
for example, for variation in metallicity between the bar and the 
disk, to do a comparison of HII region properties, and to compare 
the old population of the disk with the pseudobulge. Adding more 
stellar indicators, like the Call triplet or the NaD, would be a good 
test for the intermediate (between 10 and 100 Myr) stellar popula- 
tions. A detailed dynamical analysis of the galaxy, to directly fol- 
low the material in the disk and along the bar would be important. 
This study is also being repeated for more galaxies with different 
proprieties, Cantin et al. in preparation. 



7 CONCLUSIONS 

The integrated spectrum covering the OASIS field of view on the 
galaxy central 780pcx680pc shows strong nebular emission as 
seen in H II regions as well as absorption lines of Mg2 and Fe I 
possibly typical of an old stellar population. From this, we have 
performed an analysis using synthesis codes and an iterative tech- 
nique to study, in a consistent way, two different stellar populations 
in all the individual spectra obtained. This technique takes into ac- 
count the contribution to the flux from both populations and uses 
Bayesian statistics to get reliable uncertainties on the large number 
of parameters and spectra. We conclude that: 

• The instrumental high spatial resolution allows us to separate 
and identify the position of the young star-forming populations. 
They are distributed in a few knots on each side of the galaxy center 
and oriented along the galaxy large scale bar. A dusty nuclear bar 
is seen perpendicularly to the galaxy large scale bar. The old stellar 
population covers in a quasi uniform way the whole field of view. 

• The young bursts of star formation have an age between 5.5 
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and 8Myr with a metallicity near 2Zq. They represents nearly 40 
per cent of the stellar mass in the OASIS field of view. The old 
population age is distributed between 100 Myr and 3Gyr with a 
peak around 1 Gyr. Its metallicity is 0.2 Zq on average, a factor of 
10 lower than the gas metallicity. 

• We speculate about a simple scenario to explain this activity 
which involve slow secular evolution to build a pseudobulge for the 
galaxy. In this scenario the galaxy large scale bar, as well as the star 
forming knots and the dusty nuclear bar are playing a major role. 

• An important advantage of the iteration technique is to de- 
crease the uncertainty in the determination of age and mass of the 
stellar populations. There is a clear contribution to the continuum 
due to the younger population which needs to be taken into con- 
sideration, as we did in the iteration technique. It is also impor- 
tant to account for the spatial distribution into multiple knots of 
the younger population. For NGC 4900, the age of the old stellar 
population obtained after the final iteration is larger than the one 
estimated from iteration 0. On the contrary, the young stellar pop- 
ulation is at the end younger than obtained from the first analysis. 

• The iteration technique has limits. It assumed a mixture of 
only two populations. To some extent, if the old stellar population 
is too small relative to the young one, it becomes impossible to 
study the absorption lines and characterize the old stellar popula- 
tion. An important non-thermal contribution to the spectrum would 
not be treated properly with the actual technique. The technique 
could also be improved by considering a more detailed treatment 
of the extinction and more models for the absorption lines of H/3, 
Mg2, and Fe I in young and metal rich stellar populations. 

• Finally, the high spatial and spectral resolution of OASIS com- 
bined to spectral sjmthesis codes and our iteration technique are 
extremely powerful to extract information about the morphology, 
content, and history of galaxies. 
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APPENDIX A: BAYESIAN STATISTICS 

For the detailed analysis of § 5, we used a probability program 
based on Bayesian statistics dGregorvlliool) to extract: 1) a single 
age-metallicity solution for a stellar population when comparing its 
observational set of indicators (absorption or emission lines) with 
models, and 2) the individual mass of two populations when com- 
paring the observational spectrum with synthetic spectra. One main 
advantage of this program is to quickly give reliable uncertainties 
for a large number of spectra. We give here a brief description of 
the probability program used. 

The first step is to define the prior information W{pi,p2\I) 
based on the previous iteration I. The parameters, pi and p2, are 
either: 1) the age and the metallicity of a population, or 2) the mass 
of the old population and the mass of the young population. We 
then compute the likelihood jC{D\pi,p2, 1) of the models, which 
is the probability of obtaining the data D given that pi, p2 and / 
are true. Each model gives many indicators Li[pi,p2], which are 
related to measurements of the lines or the continuum level. The 
likelihood is calculated for TV indicators using the following equa- 
tion: 



C{D\pi,p2,I) = (27r)- 



exp 



,(A1) 



where ai is the uncertainty measured for the indicator Oi and a is 
the product of the uncertainty for all indicators. 

A posterior probability P[pi,p2\D, I] is then calculated by 
applying the prior W[pi,p2|/] to each likelihood found for all the 
models considered and by normalizing according to the equation: 



P[pi,P2\D,I] 



W[pi,p2\I]*£.[D\pi,p2,I] 



(A2) 



j:i,„,,]W[PuP2\I]*C[D\pup2,I] 
The highest probability, i.e. closest to 1, indicates the best so 
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Figure Al. The marginalized probability distribution for the age of the stel- 
lar population in the spectrum #160 in the region LI. These plots are the 
results of the Bayesian analysis: (a) for the young stellar population, and 
(b) for the old stellar population. 



lution. To compute the incertainties we marginalize the posterior 
probability over one parameter according to equation: 



p[p^\d,i] = J2p[p^,P2\d,i]. 



(A3) 



When searching for the age-metallicity pair solution for a young 
stellar population, we use a semi-flat prior for the first iteration. In 
this case, all ages are given the same weight but the models with 
the metallicity closer to the one measured for the gas is given twice 
as much weight as the other three possibilities (there are four pos- 
sible metallicities with LAVALSB). Subsequent iterations will adopt 
the posterior probability of the previous iteration as the new prior. 
Fig. Al shows, for example, the marginalized probability plot ob- 
tained for the age of the young stellar population in one spectrum 
(#160) in the region LI after the last iteration (as described in § 5, 
the iteration technique takes into account the underlying older stel- 
lar population). It was obtained by comparing two indicators, the 
equivalent widths of H/3 and Hq, with those produced by LAVALSB. 
This plot shows the probability marginalized over the metallicities, 
i.e. the probability at a given age summed over all four possible 
metallicities. The marginalized probability is used here to visualize 
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Figure A2. The probability density plot for the mass of the young and 
old stellar population in the spectrum #160 in the region LI based on the 
Bayesian analysis. 




(b) 



the solution and calculate an uncertainty for the best model param- 
eters. Vertical lines on this plot show the region containing 68 per 
cent of the total probability distribution, i.e the 1 sigma deviation 
for the best parameter pair. We see here a narrow Gaussian indicat- 
ing an age of 7.3±0.2 Myr for the young population. 

When searching for the age-metallicity pair solution for an 
old stellar population, the first iteration considers a flat prior, which 
gives the same weight to all models (there are more indicators in 
this case). Subsequent iterations are computed using the posterior 
probability found in the previous iteration. Three to four indica- 
tors are available in the spectrum for the old stellar population, the 
Fe I A 5270, Fe I A 5335, Mg2 indice, and the H/3 line when not con- 
taminated by emission (i.e after the first iteration). Fig. Al shows, 
for example, the marginalized probability plot obtained for the age 
of the old stellar population in the same spectrum (#160). We obtain 
log(age [yr])=8.4±0.5. 

When the method is used to find the mass of the young and the 
old stellar population, a flat prior is considered for the first iteration. 
It seems more efficient to keep a flat prior for the subsequent iter- 
ations. In this case, the MRl and the MR2 spectral range are con- 
sidered as two indicators (i.e. the sum of all Li — Li[young, old] 
is first calculated for each pixel of MRl and then MR2 and treated 
as two indicators in eq. Al). Fig. A2 shows the probability density 
plot obtained after the last iteration of the same spectrum presented 
in Fig. Al. Such a plot is possible since the mass parameter space 
is much larger than the metallicity parameter space. Ellipses in this 
plot are guiding lines for the mass uncertainty which is again deter- 
mined using the marginalized probability as shown in Fig. A3 for 
the young and old stellar population, respectively. 



1.0x106 1.5x106 

Mass M© 



Figure A3. The marginalized probability distribution for the mass of the 
stellar population in the spectrum #160 in the region LI . These plots are the 
results of the Bayesian analysis: (a) for the young stellar population, and (b) 
for the old stellar population. 



